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Abstract 

In this paper we use a zero-range potential (ZRP) method to model positron interac- 
tion with molecules. This allows us to investigate the effect of molecular vibrations 
on positron-molecule annihilation using the van der Waals dimer Ki2 as an example. 
We also use the ZRP to explore positron binding to polyatomics and examine the 
dependence of the binding energy on the size of the molecule for alkanes. We find 
that a second bound state appears for a molecule with ten carbons, similar to recent 
experimental evidence for such a state emerging in alkanes with twelve carbons. 
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1 Introduction 



The main aim of this work is to achieve better understanding of large anni- 
hilation rates observed for many polyatomic molecules [1,2,3,4]. In particular, 
we use an exactly solvable model to verify the prediction [5,6] that positron 
capture into vibrational Feshbach resonances (VFR) gives rise to a strong 
enhancement of the annihilation rate. We also use this model to investigate 
the dependence of positron binding energy for polyatomics on the size of the 
molecule. Such binding was postulated in [5,6] as a necessary condition for 
VFRs. 
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The annihilation rate for positrons in a gas of atoms or molecules can be 
expressed in terms of an effective number of electrons, Z ef[ , by 

A = nrlcZ cS n, (1) 

where r is the classical electron radius, c is the speed of light and n is the 
number density of the gas. Measurements by Paul and Saint-Pierre in the 
early sixties [1] indicated unusually large positron annihilation rates in cer- 
tain polyatomic molecular gases, with Z eff exceeding the actual number of 
electrons by orders of magnitude. They speculated that this might be caused 
by the formation of positron-molecule bound states, and later Smith and Paul 
[7] explored the possibility that the large rates might be caused by a vibra- 
tional resonance. Research on the alkanes and similar molecules since that 
time [2,3,4] uncovered a rapid growth of Z eS with the size of the molecule and 
very strong chemical sensitivity of Z cS . However, only recently a verisimilar 
physical picture of this phenomenon has begun to emerge [5,6,8]. These pa- 
pers put forward a mechanism which is operational for molecules with positive 
positron affinities, and which involves capture of positrons into VFRs. 

Recent measurements of annihilation with a positron beam at high resolu- 
tion (25 meV) [9,10], have shown resonances in the energy dependence of the 
annihilation rate parameter, Z e g, of alkane molecules. Most of the features ob- 
served have been unambiguously identified as related to molecular vibrations. 
In particular, for all alkanes heavier than methane Z c g displays a prominent 
C-H stretch vibrational peak. The experiments found that the magnitude of 
Z e Q in the peak increases rapidly with the size of the molecule (similarly to 
the increase in Z c g observed with thermal room-temperature positrons [2,3,4]). 
Another remarkable finding concerns the position of the C-H peak. While for 
ethane its energy is close to the mode energy of C-H stretch vibrations (0.37 
eV), for heavier alkanes the resonances appear at an energy ~ 20 meV lower 
for each carbon added. This downward shift provides evidence of positron 
binding to molecules. The binding energies observed increase from about 14 
meV in C^Hs to 205 meV in Ci2H 2 6- Very recent experiments show evidence 
of a second bound state for alkanes with 12 and 14 carbons [11]. 

So far, realistic molecular calculations have not been able to reproduce en- 
hanced Z e g. For hydrocarbons and a number of other polyatomics, calcula- 
tions have been done using a model positron-molecule correlation potential 
in a fixed nuclei approximation [12,13]. Such calculations often provide a rea- 
sonable description of low-energy positron-molecule scattering. However, their 
results, almost without exception, underestimate experimental Z eS , in some 
cases by an order of magnitude. This suggests that to describe enhanced Z e g, 
dynamical coupling to molecular vibrations must be included. Such coupling 
was considered earlier for diatomics and CO2 [14,15], where it had a relatively 
small effect on Z cS . (These molecules most likely do not form bound states 
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with the positron, and do not possess VFR.) Calculations by the Schwinger 
multichannel method [16] , which treats electron-positron correlations ab initio 
for fixed nuclei, also underestimate Z ef[ for molecules such as C 2 H 4 [17] and 
C 2 H 2 [18] by an order of magnitude [19]. 

To examine the effect of vibrations on positron scattering and annihilation, 
we consider a simple model of Kr 2 dimer using the zero-range potential (ZRP) 
method [20]. In this model the interaction of the positron with each of the 
atoms is parametrised using the atomic value of the scattering length. It is 
applicable at low energies when the de Broglie wavelength of the projectile 
is much larger than the typical size of the scatterers. Once ZRP is adopted, 
the problem of the positron-molecule interaction, including the vibrational 
dynamics, can be solved practically exactly. In the previous paper [21] the 
interaction between the atoms in the dimer was treated using the harmonic 
approximation (HA), which allowed the vibrational coupling matrix elements 
to be calculated analytically. A parabolic potential does not describe well the 
shallow asymmetric interatomic potential for a weakly bound van der Waals 
molecule such as Kr 2 . In this work we use the Morse potential to provide a 
better description of the molecular interaction. It is a good approximation 
to the best potential available for Kr 2 [22]. We examine how the use of a 
realistic molecular potential affects the positions and magnitudes of the VFR. 
To explore positron binding to polyatomics we again use the ZRP method. 
Specifically, we model alkanes by representing the CH 2 and CH 3 groups by 
ZRPs. We investigate the dependence of the binding energy on the number 
of monomers and find that a second bound state emerges for a molecule with 
ten carbons. 



2 Zero-range model for a molecular dimer 



In a van der Waals molecule the atoms are far apart and are only weakly 
perturbed by each other. This makes it an ideal system for applying the ZRP 
method. In this work we model the interaction between the two Kr atoms 
using the Morse potential (MP), 

U(R) = U miQ [e- 2a(R -^ - 2e- a ( R - Ro \ (2) 

with the parameters Ro = 7.56 a.u., C/ m i n = 6.32 x 10~ 4 a.u. = 17.2 meV, and 
uj = (2U min a 2 /m) 1 ' 2 = 1.1 x 10 -4 a.u. = 2.99 meV [23,24], where m is the 
reduced mass of Kr 2 . The energy eigenvalues and eigenfunctions of the MP 
are given by simple analytical formulae [25]. In Fig. 1 we compare the Morse 



3 



potential to an accurate fit of the best available Kr 2 potential [22], 



V(R) = Ae- aR -? R2 -J2 hn(R,b) 



a 



2d 



R 2r - 



(3) 



where a, j3, and A characterise the short-range part of the potential, and C 2n 
is a set of six dispersion coefficients. The function f 2n (R,b) is the damping 
function [26], 

{bR) k 



■in 



f 2n (R,b) = l-e- bR Y: 



k=0 



k\ 



(4) 



The values of the parameters given in [22] are: a = 1.43158, (3 = 0.031743, 
A = 264.552, b = 2.80385, C 6 = 144.979, C$ = 3212.89, C 10 = 92633.0, 
C l2 = 3.57245 x 10 6 , C 14 = 1.79665 x 10 8 , and C 16 = 1.14709 x 10 10 (atomic 
units are used throughout). 
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Fig. 1. Comparison of the best Kr 2 potential (solid curve) with the Morse potential 
(dashed curve) and harmonic approximation (dotted curve). Chain curve is the 
adiabatic potential for the e + Ki2- 

Figure 1 shows that the Morse potential is close to the best Kr 2 potential, 
while the HA is valid only in the vicinity of the minimum. This conclusion 
is supported by the comparison of the vibrational spacings. For the MP we 
have luxo = E\ — Eq — 2.74 meV, uj 2 i = 2.47 meV, uj^ 2 = 2.21 meV, which 
agree well with uio = 2.65 meV, u 2 i = 2.39 meV, 0^32 = 2.12 meV for the Kr 2 
potential [22]. Both potentials are strongly anharmonic, with u n+ i^ n deviating 
markedly from u = 2.99 meV of HA. Obviously, the MP is a much better 
approximation than HA for modelling the Kr 2 potential. 

In the ZRP model the interaction between a positron and an atom is expressed 
as a boundary condition for the positron wavefunction 



1 d(rip) 
rip dr 



-«o, 



(5) 
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where kq is the reciprocal of the scattering length [20] . Positron- Kr scattering 
calculations yield kq = —0.1 a.u. [27,28,29]. When applied to a two-centre 
problem, this condition can be expressed as 



| r ^ Ri ~ const x — - k J , (6) 

where r is the positron coordinate, and Rj (i — 1, 2) are the coordinates of 
the two atoms. 

Outside the (zero) range of action of the atomic potentials, the positron- 
dimer wavefunction can be written as a linear combination of the incident and 
scattered waves, 

„jfc„|r-Ri| „ifc n |r-R 2 | 

= e ik °' r <t> (R) + EA,$n(R), 5-r + £fl„$„(R)-j iTT> ( 7 ) 

n |r — ±ti| n |r — ±t 2 | 

where k is the incident positron momentum, $ n is the nth vibrational state 
of the molecule (n — 0, 1, . . . ), and R = Ri — R 2 is the interatomic distance. 
Equation (7) is written for the case when the positron collides with a ground- 
state molecule. The coefficients A n and B n determine the excitation amplitude 
of the nth vibrational state of the molecule, and k n = [kg — 2(E n — Eq)] 1 / 2 is 
the corresponding positron momentum. 

Applying (6) gives a set of linear equations for A n and B n , 

B m = -(e lk »^ 2 ) n0 , (8) 

/ nm 

(«o + iK)B n + -p~ ) A m = -(e- tk0 - nR/2 U, (9) 



nm 



where n is a unit vector along the molecular axis (whose direction we assume 
to be fixed during the collision), and the matrix elements are given by 

(e ±iko nR/2 ) nQ = J K(R)e ±lk °' nR/2 MR)dR, (10) 

-jp) = J K(H)—<s> m (R)dR. (li) 

/ nm 

In HA these matrix elements can be evaluated analytically, (10) - exactly, and 
(11) in the leading order [21]. For the MP we calculated them numerically. 

After solving equations (8)-(9) for A n and B n , one finds the total elastic 
(0 — > 0) and vibrational excitation (0 — > n, n — 1, 2, . . . ) cross sections, 

a ^ n = ^\A n + B n \ 2 , (12) 

Kq 
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and the positron annihilation rate, 

Z cS = ziSKlY / (\A n \ 2 +\B n \ 2 ), (13) 

n 

where zj^ is the positron-atom Z e g at k = (see [21] for details). For Kr we 
use Zgjj = 81.6 [27]. Equations (8)-(9) also allow one to determine the energies 
of bound states of the positron-dimer system, by looking for the poles of A n 
and B n at negative projectile energies, i.e., for imaginary positron momenta 
ko = i\ko\. 

For doing numerical calculations, the set of equations (8)-(9) can be truncated 
by assuming that A n = B n = for n > N c . This means that only the first 
N c + 1 channels with n — 0, 1, . . . , N c are taken into consideration. At low 
projectile energies only a small number of channels are open, and one obtains 
converged results with a relatively small N c . In the calculations we used N c = 
15 and 10, for the HA and MP, respectively. This value for MP is the total 
number of bound excited states. 

In the single-channel approximation, N c = 0, the HA results practically coin- 
cide with those of the fixed-nuclei approximation, since the matrix elements 
(10) and (11) become e ±lk »- nR o/ 2 anc i e lkmR ° / Rq, respectively (neglecting the 
2nd-order and higher corrections in the small parameter ko(muoY l l 2 [21]). A 
similar calculation for MP produces slightly different results, because of the 
asymmetry of the vibrational ground-state wavefunction, which gives rise to 
first-order corrections to these matrix elements. 



3 Results for Kr 2 

Table 1 shows the values of the bound states (negative) and the VFRs (pos- 
itive) of the e + Kr2 complex obtained with MP and in HA. In the N c = 
approximation the binding energies are eo = —3.77 meV and e = —3.48 meV 
for HA and MP, respectively. The binding energy for the MP is smaller due 
to the asymmetry of the potential curve. The corresponding energies obtained 
from a multichannel calculation given in Table 1 are lower, because allowing 
the nuclei to move leads to stabilisation of the e + Kr 2 complex. 

The ground-state energy of the complex can also be compared to the results 
of an adiabatic calculation. For fixed nuclei the energy of the positron bound 
state is —k 2 /2, where k is a positive root of the equation k — kq + e~ KR /R. 
Adding this energy to the Kr 2 potential, one obtains the adiabatic potential 
for the e + Kr 2 complex (chain curve in Fig. 1). Its minimum is about 3.94 meV 
below that of the Kr 2 , which is close to the MP value of e in table 1. 
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Table 1 

Energies of the bound states and resonances for e + Kr2. 





Energy 


(meV) 


Level 


HA 


MP 





-4.23 


-3.89 


1 


-1.41 


-0.74 


2 


1.42 


2.16 


3 


4.25 


4.83 



The first vibrational excitation energy of the e + Kr2 complex, uj[ = E\ — Eq, for 
MP is 3.15 meV, while in HA it is 2.82 meV. Thus, MP calculations predict a 
"stiffening" of the vibrational motion of the complex in comparison with that 
of Kr 2 . This effect is caused by a shift of the equilibrium position of the atoms 
to the left (see Fig. 1), towards the steep repulsive part of the interatomic 
potential. Note that in MP, the 2nd bound state with E\ = —0.74 meV lies 
just below the threshold. We will see that this causes a steep rise in o"o-+o an d 
Z e g as k — > 0. In HA this bound state is lower, i.e., further away from threshold, 
and its effect on the cross sections is less noticeable. This combination of a 
lower binding energy and a greater vibrational frequency in MP, means that 
the first resonance observed in the cross sections and in Z e s will be at a greater 
energy than in HA. 

Figure 2 shows the elastic and vibrational excitation cross sections obtained 
from Eq. (12) for MP and in HA. To highlight the effect of resonances, the 
elastic scattering cross section has been calculated in both multichannel and 
single-channel (N c = 0) approximations. The single-channel ("fixed nuclei") 
elastic cross sections from the two calculations are quite close. The multichan- 
nel cross sections are qualitatively similar, the main difference being in the 
positions and widths of the resonances and energies of the vibrational thresh- 
olds. The magnitudes of the a ^i and <t ^2 cross sections are greater for MP. 
Another noticeable difference is in the rise of towards zero positron en- 
ergy in MP calculation. 

Figure 3 shows the positron annihilation rate (13) obtained with and without 
the coupling to the vibrational motion, i.e., from the multichannel and single- 
channel calculations. The background ("fixed nuclei", N c = 0) Z cS at low 
positron momenta is enhanced due to the large positron-Kr 2 scattering length. 
Such enhancement first predicted in Ref. [30], affects both Z cS and the elastic 
cross section, which in this case are proportional to each other, Z e g ~ o" e i/47r 
in atomic units [5,6]. The effect of VFRs on Z e g is much more prominent than 
in scattering, with the strongest resonance four orders of magnitude above the 
background. The widths of the resonances in MP and HA, are quite different, 
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Fig. 2. Cross sections for positron scattering from Kr2 calculated using HA (left) 
and MP (right): solid curve, elastic scattering, cn^oi long-dashed curve, vibrational 
excitation, <7n— >i (times 10 2 ); chain curve, vibrational excitation, <7n— >2 (times 10 4 ); 
short-dashed curve, single-channel elastic cross section. 

e.g., T = 2.8 fieV (MP) vs. T = 16.7 /j,eV (HA) for the strongest n = 2 
resonance. This difference, also seen in the scattering cross sections (Fig. 2), 
means that anharmonicity of the Kr2 potential reduces the coupling between 
the incident positron and vibrationally-excited e + Kr 2 compound. A possible 
explanation for this is that positron binding has only a small effect on the 
equilibrium position of the nuclei (as seen from adiabatic potential curve of 
e + Kr2 in Fig. 1). 
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Fig. 3. Z e Q for positrons on Kr2 calculated using MP and HA in the multichannel 
(solid and chain curves, respectively) and single-channel (dashed and dotted curves, 
respectively) approximations. 

To compare the integral contribution of the resonances, we averaged Z e ^ over 
the Maxwellian positron energy distribution, 



Z eS (T) 



e -k 2 /2k B T 



(14) 



Figure 4 shows Z e g(T) for HA and MP. In both cases the VFR gives a very 
large contribution, increasing Z e g by an order of magnitude at T ~ 1 meV, 
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i.e., for positron energies close to that of the resonance. Its effect is seen even 
at much higher temperatures, raising Z e g above the non-resonant background 
by a factor of three for room temperature positrons. 




Temperature (meV) Temperature (meV) 



Fig. 4. Thermally averaged Z e g for positrons on Kr2, obtained using HA (left) and 
MP (right): long-dashed curve, single-channel approximation; chain curve, multi- 
channel calculation, non-resonant background; solid curve, multichannel calcula- 
tion, total, including the VFR. For comparison, Z c r for Kr is shown (short-dashed 
curve) . 



4 ZRP model of positron binding to polyatomics 

In positron-molecule collisions VFRs occur when the energy of the incident 
positron plus the energy released in positron-molecule binding, equals the en- 
ergy of a vibrational excitation of the positron-molecule complex. For weakly 
bound positron states the latter should be close to a vibrational excitation 
energy of the neutral molecule. Hence, by observing VFRs one can obtain 
information on the binding energy. This procedure was applied to electron 
scattering from (CC^n clusters [31]. In this system the redshift of the VFR 
was found to increase with the cluster size by about 12 meV per unit. A sim- 
ple model of a spherically-symmetric cluster with a constant potential inside 
and — oe 2 /r 4 potential outside was able to reproduce the dependence of the 
electron binding energy on the cluster size. 

In a similar way, the measurements of the energy dependence of Z e ff for alka- 
nes with a high-resolution positron beam allow one to determine their positron 
binding energies [9,10]. In contrast, an accurate ab initio calculation of positron 
binding to a polyatomic molecule is probably beyond the capability of present- 
day theory. Even for atoms, calculations of positron binding remain a challeng- 
ing task because of the need to carefully account for strong electron-positron 
and many-electron correlations (see, e.g., [32]). 

In this work we have adopted a different approach. To examine positron bind- 
ing to alkanes, we model the molecule by representing each CH 2 or CH 3 group 
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by a ZRP centred on the corresponding carbon atom. The wave function of 
a bound state decreases as r _1 e~ Kr at large positron-molecule separation r, 
where k 2 /2 is the binding energy. For weakly bound states this wavefunction 
is very diffuse (k < 1 a.u.), which means that the positron moves mostly 
far away from the molecule. The actual binding energy is determined by 
their interaction at small distances, and the ZRP model is a simple way of 
parametrizing such interaction. It allows us to account for the scaling of the 
positron-molecule interaction with the number of monomers (to the extent 
that the monomers do not have a large effect on each other), and to use a 
realistic geometry of the molecule. We will consider two straight car- 

bon chain, Fig. 5 (a), and a "zigzag" carbon chain, Fig. 5 (b), which matches 
the actual structure of the molecule, Fig. 5 (c). Unlike the Kr 2 model, the k 
parameter of the ZRP for alkanes is adjusted semiempirically (see below). 




(a) 

(b) 

00 

Fig. 5. Zero-range potential models of the alkane molecule (butane, C4H10, is 
shown). In (a) and (b) the shaded circles represent a CH2 or CH3 group, while 
(c) is a true 3D molecular structure [33]. The parameters used are Rq = 2.911 a.u. 
and 9 = 113°. 

The bound-state positron wavefunction in the field of iV ZRP centres has the 
form [20], 

N „-k|t— Ri| 

* = £At— ^j. (15) 
i=i 11 rt »- 



Subjecting it to iV boundary conditions (6) with parameters K 0i , we find the 
positron energy as Eq = — k 2 /2, where k is a positive root of the equation 



det 



(K 0i - K)5ij H — - — (1 - 5ij) 



0. (16) 



Here = |Rj — Rj| is the distance between the ith and jth ZRP. 

For modelling alkanes we choose the distance between the neighbouring ZRPs 
equal to the length of the C-C bond. All ZRPs are characterised by the same 
value of K,Qi = —1.5. This value ensures that the molecule with three ZRP 
centres (which models propane, CsHg) has a small binding energy, 7 meV for 
the straight chain, and 12 meV for the zigzag chain. These values are close to 
that inferred from experiment [9,10], where propane is the first molecule for 
which a downshift of the C-H peak from the corresponding vibrational energy 
can be seen [9,10]. 
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In Fig. 6 the results of our calculations are compared with the experimental 
binding energies. As we move from a straight ZRP chain, Fig. 5 (a), to a 
"zigzag" chain, Fig. 5 (b), the binding energy increases. This is expected as 
the carbon atoms beyond the nearest neighbour become closer to each another. 
The overall dependence of the binding energy on the number of monomers n 
predicted by the ZRP model is similar to that of the experimental data, while 
the absolute values predicted by our simple theory are within a factor of two 
of the measurements. One may also notice that the measured binding energies 
increase almost linearly with n, while the calculation shows signs of saturation. 
These discrepancies might be related to the fact that the ZRP model disregards 
the long-range — oe 2 /r 4 nature of the positron-molecule interaction, which 
would restricts its validity to very small binding energies. 
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Fig. 6. Positron binding energies |eo| for alkanes modelled using straight (circles) 
and "zigzag" (triangles) ZRP chains. Experimental results (crosses) [9,10] are shown 
for comparison. 



A remarkable feature of the model calculations is the emergence of a second 
bound state for n = 10 in both straight and "zigzag" chains. This prediction 
supports the experimental evidence for the second bound state, in the form 
of a C-H peak, which re-emerges at 0.37 eV for dodecane (n = 12) [9,10] and 
stabilises by about 50 meV for C14H30 [11]. 



5 Summary and outlook 



We have used the Morse potential to model the interaction between the atoms 
in the Kr 2 dimer. We find that the positron binding energies and the positions 
and widths of the VFR change compared with the harmonic approximation. 
However, the overall picture remains similar, with the lowest VFR enhancing 
the Maxwellian-averaged positron annihilation rate by an order of magnitude 
to Z cS ~ 10 4 at T ~ 1 meV. 
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In priniciple, a similar approach could be applied to positron interaction with 
other rare-gas dimers and clusters. For Ar and lighter atoms, the positron- 
atom attraction is too weak to allow formation of positron bound states and 
VFRs. For Xe 2 , on the contrary, the attraction is much stronger (/c ~ —0.01 
[28]). This leads to a much greater positron-dimer binding energy (~ 40 meV), 
which means that many vibrationally excited states of e + Xe2 are bound, and 
only those with high n lie above the positron-dimer continuum threshold. The 
coupling between these states and the target ground state is extremely weak, 
and we have not been able to find any VFR for positrons incident on ground- 
state Xe 2 in our calculations. 

A zero-range potential model for positron binding to alkanes yields binding 
energies that are in qualitative agreement with experiment. Our calculation 
predicts the emergence of the second bound state for a molecule with ten car- 
bon atoms. Such a bound state may have already been observed for dodecane 
and tetradecane. 

Zero-range potential is an exceptionally simple model. The accuracy of our 
predictions is of course limited by the nature of the ZRP model. In particular, 
ZRPs disregard the long-range character of the positron-target polarization 
attraction. This is a reasonable approximation for very weakly bound states, 
but as the binding gets stronger, larger errors may be introduced. Given the 
difficulty of performing ab initio calculations for positron interaction with 
polyatomics, one hopes that more sophisticated yet tractable models could 
be developed. They should provide a more accurate description of positron- 
molecule binding and capture into vibrational Feshbach resonances. One may 
then hope to fully explain the dramatic enhancement of the annihilation rates 
and their strong chemical sensitivity for polyatomic molecules. 
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